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Using an inductive approach to classify multipartite entangled states under stochastic local op- 
erations and classical communication introduced recently by the authors [Phys. Rev. A 74, 052336 
(2006)], we give the complete classification of four-qubit entangled pure states. Apart from the ex- 
pected degenerate classes, we show that there exist eight inequivalent ways to entangle four qubits. 
In this respect, permutation symmetry is taken into account and states with a structure differing 
only by parameters inside a continuous set are considered to belong to the same class. 



PACS numbers: 03.67.Mn, 03.65.Ud, 02.10.Yn 
I. INTRODUCTION 

Entanglement is the distinguishing feature of quantum 
systems [H, 0] and is especially useful in using such sys- 
tems to process information 0, 0] ■ Despite its relevance 
this property is not yet fully understood and different 
questions remain open [f| . Among them a complete clas- 
sification of entangled states under adequately chosen cri- 
teria of equivalence stands as a formidable task. Two 
of these criteria appear as outstanding, namely the so- 
called LU and SLOCC equivalences. The former stands 
for Local Unitary, which is mathematically expressed by 
the relation |?/>) ~it/ \<t>) ^ |V0 = ® . . . \(f>) 
for unitary matrices U^ k \ whereas the latter stands for 
Stochastic Local Operations and Classical Communica- 
tion, which is mathematically translated as \ip) ^slocc 
\4>) \tp) = ® ...F^\4>) for non-singular matri- 
ces F*- k '. Although these entanglement equivalence def- 
initions can be straightforwardly generalized to include 
mixed states, we will restrict ourselves to the pure state 
case. Furthermore we will center upon the identification 
of entanglement classes under the second criterion. 

The SLOCC equivalence (also known as local filter- 
ing operations) was introduced in Q and enjoys a clear 
physical motivation: two multipartite pure states \ip) and 
\4>) are equivalent under SLOCC if one can be obtained 
from the other with non-null probability using only local 
operations upon each subsystem and classical communi- 
cation among the different parties. In this manner, the 
complete classification of entanglement for N = 3 qubits 
can be derived, namely, any genuinely entangled three- 
qubit systems is equivalent either to the GHZ state or 
to the W state. The difficulties for more general cases, 
i.e. for N > 4, was also pointed out in Ref. [f|, where 
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a non-enumerable set of entanglement classes would be 
required. 

In this work, using our preceding inductive ap- 
proach we give the complete classification of entan- 
glement for N = 4 qubits. Although this was already 
apparently performed in Ref. [1] , we argue that their enu- 
meration does not follow the same philosophy as in the 
seminal paper by Diir et al. @, In Sec. [HI we compile 
the preceding relevant results which will be later used 
in Sec. IIIII to obtain in detail all entanglement classes 
under SLOCC. Our concluding remarks are presented in 
Sec. ED 



II. REVIEW OF THE INDUCTIVE METHOD 

The forthcoming entanglement classification for N 
qubits is based on the inductive approach which the 
present authors introduced in Ref. 0] and in some auxil- 
iary results which were proved therein. The main idea of 
this approach was to establish a classification of the right 
singular subspace of the coefficient matrix of each pure 
state \ip) expressed in an adequate product basis (typi- 
cally the canonical -computational- basis). This classifi- 
cation must be carried out according to the entanglement 
classes for N —1 qubits (hence the term inductive) . Con- 
sequently, knowing in advance that there are only two 
entanglement classes for 2 qubits, namely the product 
class 00 and the entanglement class the right singu- 
lar subspace 233 of the coefficient matrix of the state of 
any three-qubit system results to be of six possible types, 
i.e. two one-dimensional subspaces 233 = span{|0-0)} and 
233 = spanjl^)} driving us to the degenerate classes 000 
and 0i\I/23, respectively, two two-dimensional product 
subspaces 233 = \<j>) ® C 2 and 233 = C 2 ® \iJj), driving 
us to the degenerate classes 02^/13 and 03^12, respec- 
tively, and, finally, another two two-dimensional classes 
233 = span{00,00} and 23J = span{00, driving us 
to the well-known GHZ and W genuine entanglement 
classes. Let us recall the reader that 233 = span{00, "J"} 
stands for the fact that the subspace 233 contains (up 
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to normalization) only one product vector (i.e. only one 
of type 00). This will sometimes be also expressed as 
2U = span{|0<£), I*)}. This convention will be profusely 
followed in the next section. We will employ these six 
tripartite entanglement classes to find those of N = 4 
qubits. 

Moreover, in the proof, we will make use of the follow- 
ing results, already proved in Q. 

Proposition 1. Let 2U be a two-dimensional subspace 
in C 2 <g> C 2 . Then 2U = span{ | (/></?) , 1*)} if, and only 
if, 20 = span{ | <f>ip) , \4>(p) + | </></?)}, where ~ denotes linear 
independence. 

Let us also recall that only two options are left for the 
structure of any two-dimensional subspace (Toj |. namely 
W = span{00,00} and 2U = span{00, where the 
above convention has been used. Complementarily, the 
following criterion to discern whether a given tripartite 
pure state belongs to any of the six entanglement classes 
will be also necessary. We previously need the following 
definition. 



Definition 1. Let 



C^=C m = ( 



Clll CH2 C121 C122 
C211 C212 C221 C222 



la) 



C (2) = c = (cm cxia Pan c 212 . (lb) 

1 V Ci21 Ci22 C221 C222 / 



C( 3 )=C 3 |12 = ( 



Cm Ci21 C2U C221 
Cll2 C122 C212 C222 



(lc) 



be the coefficient matrices of the partitions 1|23, 2 1 13 
and 3 1 1 2 , respectively, of the tripartite pure state |*) = 
Sijfe Cfjfc |* — 1} ® \j — 1} ® |fe — 1). We wiZZ define the 
matrices Wi and W2 as 



Wi = 



cm C112 
C121 C122 



W 2 = ( 0211 C212 I . (2) 
\C221 C222/ 



With these definitions and remembering that the spec- 
trum of a matrix A is denoted as o~(A), in [3] we proved 
the following 

Theorem 1. Let |$) denote the pure state of a tripartite 
system and C^ 1 ' its coefficient matrix according to the 
partition i\jk. Then 

1. |$) belongs to the 000 class if, and only if, r(CW) = 
1 for all i = 1, 2, 3. 

2. |\&) belongs to the 01^23 class if, and only if, 
r(C^) = 1 andr(C^) = 2 fork = 2,3. 

3. |\&) belongs to the 02*13 class if, and only if, 
r(C^) = 1 andr(C«) = 2 for k = 1,3. 

4- |*) belongs to the 03*12 class if, and only if, 
r(C( 3 )) = 1 andr(C^) = 2 fork = 1,2. 

5. |*) belongs to the GHZ class if, and only if, one 
of the following situations occurs: 



i. r(CW) = 2 for all i = 1,2,3 and r(Wi) = 
r(W 2 ) = 1. 

M. r(C«) = 2 for all i = 1,2,3, r(Wj) = 2, 
r(W fc ) = 1 and c^W^Wj.) «s non- degenerate. 

Hi r(CW) = 2 /or aZZ i = 1,2,3, r(Wi) = 2, 
r(W 2 ) = 2 and ^(Wf 1 ^) 

is non-degenerate. 

6. |*) belongs to the W class if and only if, one of 
the following situations occurs: 

i. r(CW) = 2 for all i = 1,2,3, r(Wj) = 2, 
r(Wfc) = 1 and a{WJ x Wk) is degenerate. 

ii. r(CW) = 2 /or all i = 1,2,3, r(Wi) = 2, 
r(W-2) = 2 and a{W^ W2) is degenerate. 



III. FOUR-QUBIT CLASSIFICATION 

The classification inductive procedure for N = 4 qubits 
is long, although systematic. We will firstly detect the de- 
generate classes, which furthermore are predictable in ad- 
vance: 0000, 0, 0,,*, , ; . O^GHZ, n W and * ili2 * i3i4 , 
with ik — 1,2,3,4. Then, we will work out the gen- 
uine classes. The classification is performed accord- 
ing to the structure of the right singular subspace of 
the coefficient matrix of the four-partite pure states: 
2U = span{* t , *.,}, where % itj = 000, fc *, GHZ, W. 



A. Degenerate classes 

As stated we will revise in turns each of the possible 
structures of the right singular subspace 2U of the coeffi- 
cient matrix. Remember that now 2U C C 2 <E> C 2 <g> C 2 . 

1. 2JJ = span{000}.- In this case, we can write 2tJ = 

span{|0<p , 0)}. Mimicking the same arguments as 
in (see also below for the genuine classes), we 
obtain the canonical vector 1 0000) , corresponding 
to the degenerate class 0000. 

2. 2JJ = span{0i*23j-.- This drives us to the canonical 

vector 1 0000) + |0011), corresponding to the degen- 
erate class 0i02 l J , 3 4. 

3. 2U = span{02^E'i3}.- This drives us to the canonical 

vector 1 0000) + |0101), corresponding to the degen- 
erate Class 0i03*24- 

4. 2U = span{ 03*12}.- This drives us to the canonical 

vector 1 0000) + |0110), corresponding to the degen- 
erate Class 0104*23- 

5. 2U = span{ GHZ}.- This drives us to the canonical 

vector 1 0000) + |0111), corresponding to the degen- 
erate class OiGHZ. 

6. 2H = span{W}.- This drives us to the canonical vec- 

tor |0001) + |0010) + |0100), corresponding to the 
degenerate class OiW. 
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7. 2U = span{000,000}.- When 2U = \(fxp} ® C 2 , we ob- 

tain the canonical vector |0000) + 1 1001) , corre- 
sponding to the degenerate class 0203*14. 

8. 2U = span{000, 000}.- When W = \4>) <g> C 2 ® |V>), we 

obtain the canonical vector |0000) + 1 1010) , corre- 
sponding to the degenerate class 0204*13. 

9. 2XJ = span{000,000}.- When 2H = C 2 ® |^>), we ob- 

tain the canonical vector |0000) + 1 1 100) , corre- 
sponding to the degenerate class 0304*3/12 . 

10. 2U = span{000,_000}.- When 22J |0) ® 

span{|</?V)> |f^V')}> we obtain the canonical vector 
|0000) + 1 1011) , corresponding to the degenerate 
class 2 GHZ. 

11. 2H = span{000 L 000}.- When 3B 

sp&n{\(jxpip) , \4>(pip}}, we obtain the canonical 
vector |0000) + 1 1 101) , corresponding to the 
degenerate class O3GHZ. 

1 2 . an = span{ 000 , 000} .- When 3D 

span{\4>(ptp) , \4><pip)} , we obtain the canonical 
vector 1 0000) + 1 1110) , corresponding to the 
degenerate class O4GHZ. 

13. 2U = span{000,0i*23}-- When 2JJ = \<f>) ® 

spaii{|(pi-0i), I "02 ) } j we obtain the canonical 
vector |0001) + |0010) + |1000), corresponding to 
the degenerate class O2W . 

14. = span{000, 2 *i 3 }.- When W 

spaii{|(/>i(^Vi), \<h l ( n h) + 102^2)}, we obtain 
the canonical vector |0001) + |0100) + 1 1000) , 
corresponding to the degenerate class O3V17. 

15. aiJ = span{000, 3 *i 2 }.- When an 

span{|^i(^i), |0 2 (p 2 ) + 102^2)} ® \ip), we ob- 
tain the canonical vector |0010) + |0100) + 1 1000) , 
corresponding to the degenerate class O4W. 

16. an = span{0i*i 3 ,0i*i 3 }.- When an = C 2 g> 

span} I 'I')}, we obtain the canonical vector 1 0000) + 
|0011) + |1100) + |1111) = (|00) + |11)) (|00) + |11)), 
corresponding to the degenerate class *i2*34- 

17. an = span}0 2 *i3, 02*13}.- When _ 2n 

span} I (jxpip) + 1 4>ifip) , 1 4"f' l P) + 1 4"f' l P) } ■ we obtain the 
canonical vector |0000) + |0101) + |1010) + 1 1111) = 



This matrix corresponds to the canonical state |0000) + 
1 1111), i.e., to the so-called GHZ state. 



(| 00) + |ll))i 3 (|00) + |11)) 24 , corresponding to the 
degenerate class * 13*24 ■ 

18. an = span}0 3 *i3,03*i 3 }.- When an = |*) ® C 2 , 
we obtain the canonical vector |0000) + 1 01 10) + 
11001) + 11111) = (|00) + |11)) 14 (|00) + |11)) 23 , 
corresponding to the degenerate class *i4*23- 

No further degenerate classes can be found within the 
rest of possible structures for an. Notice how all of them 
can be obtained from the first six and the sixteenth cases 
by just permuting the qubit indices. Wc have preferred to 
make use of the systematic approach in order to illustrate 
its usage. 

B. Genuine classes 

In order to be exhaustive we must check out 
all possible structures that the right singular sub- 
space 3n may have. An a priori list is given 
by span}000,000}, span}000, ix *}, span}000, GHZ}, 
span}000,W}, spamjO^*, l2 *}, span-fO^*, GHZ}, 
spanjo^*, 1 ^}, sp<m{GHZ,~GHZ}, spa,n{GHZ, W} 
and span}VT, W}, where %k = 1,2,3. The same conven- 
tion as above has been followed when we say that, e.g., 
spanjO^*, GHZ} indicates that only one vector belong- 
ing to the class 0^ * is contained in an, the rest being all 
GHZ (and possibly W). In other words, in span}*j, *&} 
it is understood that *^ and *& are those vectors in 30 
with the lowest degree of entanglement degeneracy, ac- 
cording to the scale 000 < tfc * < GHZ, W. 



1. 2U = span{000, 000} 

All cases leading to degenerate cases have been consid- 
ered above. However, when 3n = spaxi{\(f>ipt(} , 4>(pip)}, a 
well-known entanglement class arises, namely the GHZ 
class. We choose non-singular transformations so 
that |0) = F^\<j>), |1) = F^\4>) and similarly for 
ir[3,4] _ We choose also F^ so that | Vi ) F^ |v 2 )] = 

[mIj] -1 ( CT q 0l )' where \vj) denotes the left singular 

vectors of the coefficient matrix Cg, of the state |<J>) G 
C 2 ®C 2 <g>C 2 <g)C 2 and the coefficients expressing 

the right singular vectors |wjt) of C$ in terms of \(f>ipip) 
and \6(pip): |w;) = mi\<ptpip) + ^^(pip) . Then, we can 
write 0| 



(3) 



2. W = span}000, 0i* 23 } 

In general, in this case 3n = span{|0iy>iV>i), |(fo) ® 
(Iv^V^) + \V2^2))}- The degenerate case takes place 



J 



r (n*\- l ( ai °\ 1 ( ai 000000 fx* 12 \ 

u FMQ-9iM\*) ~ W y a2 J y (72 J ^21 M *J 

fl o\ 

^OOOOOOOlJ" 

I 
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when \4>i) and |</>2) are linearly dependent, driving us 
to the class 2 W, already stated above. 

Thus, let us consider the case when {|0i),|02)} are 
linearly independent. Now the two-dimensional sub- 
space span{|(/5i , 0i), \cp2ip2} + ^2^2)} can be either of type 
span{00, 00} or of type span{00, 

i. When span{ | }, | (p 2 i>2) + \<P2i>2)} = 



span{00,00}, then there exist linearly inde- 
pendent vectors \4>), \4>), \<ptf>) and \(p , 4>) such 
that |w fe ) = + HkiWj <8 (\<pip) + \<f>4>})- 

Thus, it is immediate to realize that there exist 
non-singular transformations F^ k \ k = 1,2,3,4, 
such that 



r („*\- 1 °\ 1 ° \ M12 

^FW 9 -9Fwm - y a2 j y U2 ) ^ o o ^ 2 o o 

(1 o\ 

^0 1 lj 



This canonical coefficient matrix corresponds to the l w fc) = Hki\4"P' l P) + A i fe2 1 0) <£> (lyV'} + W^))- 

state 1 0000) + 1 1 100) + 1 1 1 1 1 ) . Thus, it is immediate to realize that there exist 

non-singular transformations F<- k \ k = 1,2,3,4, 
ii. When spanjl^i), \<P2ih) + ^2^2)} = such that 

span{00, then there exist vectors \<f>), 

|0), \<p), \(p), \tp) and |V>) such that 



_ , r i 0\ 1 /ax 0\ //^ ^ 2 0\ 
Gf|11 ^'« p[4i|$ ' ~ [0 a 2 ) [0 (72 J ^21 fj,* 22 fj,* 22 Oj 

(l 0\ 

_ ^0 1 1 oj 



This canonical coefficient matrix corresponds to the 
state |0000) + |1101) + 1 1110) . 



3. 23 = span{000, 02*13} 

Although it is always possible to reproduce a similar 
argument to the one above, we will resort to permuta- 
tion symmetry, thus leading to the degenerate class O3W 
(already stated) and the genuine classes with canonical 
states given by |0000) + |1010) + 1 1111) and |0000) + 
|1011) + 1 1110) . 



5. 23 = span{000, GHZ} 

This class is highly richer than those above. In this 
case 2B = span{|0i<piV>i)) 102^2^2) + I^^V^)}, with 
the restriction that no product state other than l^iyiV'i) 
and no 0fc\E' state belong to 2U (otherwise we would be 
in one of the preceding classes). We define now non- 
singular transformations such that F^lfo) = |0), 
F} 2] \4>2) = |1), and similarly for \<p 2 ), \(p 2 ), H2) and 
IV2). We define also F^ such that [F^ | vi) F^ | v 2 >] = 

( 1 -1 )j with the same conventions as above. 

[ — h 

Then, the canonical vector can be written as [llj 



4. 23 = span{000, 03*12} 

With analogous arguments, in this case the degen- 
erate class is O4W (already stated) and the genuine 
classes are identified by the canonical states given by 
|0000) + |1001) + |1111) and |0000) + 1 1011) + |1101). 



|0^V) + |1000) + |1111). (6) 

The factor vectors \<f>), \tp) and \ip) are arbitrary up to 
the restriction of not producing more than one 000 vector 
and no Ofe* vector in 22J. Although this recipe embraces 
all possible cases, we will single out two different subsets 
within this class, namely (i) that of states leading to right 
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singular subspaccs 2D with no W state in it, i.e., contain- 
ing only one product vector and GHZ vectors and (ii) 
that of states leading to right singular subspaccs 2D with 
at least one W vector in it. Indeed, we will show that 
there will be either one or two W vectors in 2D. 

i. The set of states leading to 2D containing no W vec- 
tors is that identified by the canonical ones given 
by |001V>} + |1000} + |1111), |0(V1} + |1000} + |1111) 
and |0</>01) + |1000) + |1111), with no restrictions 
upon the factor vectors. Moreover, if we denote the 
components of \(f>) in the (non- necessarily orthonor- 
mal) basis {|02) 5 102)} by 0o and 0i and similarly 
for \(p) and then whenever 7^ \/0oW^o = 
iV^iViV^i 7^ 0, there will also be a unique product 
vector, the rest being all GHZ in the right singular 
vector of the generic state ([6]) . 

This is proved resorting to thcorcmQ] We will illus- 
trate the calculation with the first canonical vector 
|001^) + 1 1000) + |11H), whose right singular sub- 
space is spanned by |O10) and ^-|000) + -^=|111). 
Thus, a generic vector in 2D will be of the form 
a(|000) + |111)) + /3|O10), whose coefficient matrix 

is (cf. dD-cg)) 



(1) _/ a o« m 

lo 



a 



(7) 



For any values of tpl it is clear that r(C^) = 2 for 
« = 1,2,3 (except for the trivial case a = 0). Then 
if Vi = 0, r(Wi) = r(W 2 ) = 1, thus it is a GHZ 
vector; if Vi + 0, r(Wi) = 2 and t(W 2 ) = 1 and 
the spectrum of W^ 1 W 2 is never degenerate. It is 
again a GHZ vector for any value of \tj)). 

The other canonical vectors are handled in a sim- 
ilar fashion and they can additionally be obtained 
under permutations among qubits 2, 3 and 4. 
Among them we single out the particular cases 
|0000) + 1 1101) + |0111), |0000) + |1011) + |0111) 
and |0000) + 1 1110) + |0111), which can also be ob- 
tained from those in the O^'l' classes by permuta- 
tions among the four qubits. 

For the generic state ([6]), let us write, as explicitly 
supposed, 2D = span{|0<^>), [000) + [111)}, with 
no vector |0), and \tp) equal to |0) or |1). If we 
denote again by subindices and 1 the components 
of each vector |0), \tp) and \ip) in the {|0), |1)} basis, 
then the coefficient matrix of a generic vector in 2D 
will be 



where no coefficient is null, as supposed. Then 
it should be clear that r(C^) = 2 for all i = 
1,2,3. Furthermore, since detWi = af34>o^p\tp\ 
and detWb = ctf3(f)iip ip , we will have r(FFi) = 
T {W 2 ) = 2 (except for the trivial cases a ■ j3 = 0), 
thus we must check the degeneracy of the spec- 
trum of W^ X W 2 : a{W^ x W-i) is degenerate when- 
ever the discriminant of the characteristic poly- 
nomial is null, i.e. whenever (tr(W 1 -1 W2)) — 
4det(W 1 ~ 1 W , 2) = 0, which drives us to the con- 
dition f3 + a (v / 0o¥>o'0o ± V^hv^h) 2 = 0, which 
immediately yields the above stated conditions. 

ii. The set of states leading to 2D containing at least 
one W vector comprises the rest of possibilites. To 
prove this assertion we need to show that all pos- 
sibilities for 2D having no W vector has been al- 
ready considered. The most generic case occurs 
when 2D = span{|0i<£iV>i), 102^202) + 102^202)}, 
which we detach in the following particular cases: 

• If \<j)i<piipi) = I02V2020 10202020 then two 
000 vectors belong to 2D. This case is ruled 
out. 

• If Ifatpilpl) = |0 2 V?20), 102020), With |0) jt 

1^2)) 1020 respectively, then there will be 
a Ofc^ in 2D unless \tp) = $2), 1020 re- 
spectively. We give the proof for the first 
case. Let |0iv?i0i) = |02<P20), with |0) ^ 
|020 A generic vector in 2D will be of the 

form a\<f>2<P2lp) + P (|02<£202) +_ |02<^2'02>) = 

102^2) H0) + /3\tp 2 )) + 0|0202_020 It is clear 
that, provided that |0) and \ip 2 ) are linearly 
independent, it is always possible to find a, /? 
such that a\tp) + P\ip 2 ) = P\ l p2), driving us 
to a 03^ vector in 2D, against the general hy- 
pothesis. By symmetry, all other cases can 
be accounted for in a similar fashion, leaving 
us in the cases that whenever two vectors in 
\<j>tpip) _are equal to two vectors in \<p 2 ip 2 ip 2 ) or 
10202020 then the other must be the third one 
corresponding to the opposite. For instance, 
if 102^202) = |000), 102^2-02) = |111) and 
Ifatpiipx) = |OO-0), it must be \ip) = |1). 

• If |0i<^i0i) = |0 2 ^0),_with \ip) ^ \ip 2 ), \ip) ^ 
\tp 2 ) and \<pip) 7^ 102020 use the same non- 
singular F*- k ' as in the generic case (|6]) to arrive 
at the coefficient matrix for a general vector 
in 2D 



c*vo0o + P a(p ipi aip^o oifiipi 
/3 



(9) 



C 



' a0 o <^o0o + (3 afotpoip! a0o<Pi0o a0o<?i0i 
^ a^i^oV'o a0i<^ o 0i afatpx-ipo a0i<£i0i + 

(8) 



with the above conditions translated as \tp) 7^ 
|0) «^i^0, |V) |0) & 0i and \tp^) ± 
|11) (tpo,ifjo 7^ 0) simultaneously. It is then 
clear that (except for the trivial cases a-(3 = 0) 
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r(C«) = 2 for all i = 1,2,3, r(Wi) = 2 and 
r(W 2 ) = 1. The spectrum a{W^ 1 W l ) is de- 
generate only if aifiQipo +(3 = 0. Thus there 
will be no W vector in 20 only when foipo = 0, 
i.e. when \cp) = |1) or \tjj) = |1), which going 
back before the application of the non-singular 
FV°>, means that \tpi) = \(p2) and |V>i) = l^)- 
By using <r x upon qubits 2,3,4 the symmet- 
ric case_|</>i^i"0i) = \(fo<pi>), with \<p) ^ \tp 2 ), 
H') 7^ I "02) and \(ptp) ^ \ f2ip2) is also included 
in this analysis. By permutations among the 
qubits 2, 3, 4 all preceding cases are accounted 
for in a similar fashion. 

• Finally the case l^i^'iV'i) — \4>ip%j)), with 
\4>) ±]<h),\h)\ \<p) ± |^2>, |^2> and M ± 
^2) ,^2) , drives us, after application of the 
non-singular , to the coefficient matrix (|8|), 
already accounted for. 

No more options are left, thus we have scrutinized 
all possible 2U's with structure span{000, GHZ}. 



6. 22T = span{000, W} 

By span{000, W} we indicate that only one 000 vector 
and W vectors belong to 2U. The generic case will be that 
expressed by 2U = span{|0i^i'i/'i}, 102^202) + |</>2 ^2^2} + 
l&^fa)} ■ We will show that there will be no GHZ state 
in 2U only if \<piLpiipi) = ^2) , which after applica- 
tion of the non-singular transformations F^ will drive us 
to the canonical state |0001) + |0010) + |0100) + |1000), 
i.e. the W state for 4 qubits. As before, we will be ex- 
haustive: 

• If = |</>2</?2?/>2), 102^2^2)) 102^2^2), then 

it is clear that a Ofc^ can be found in 22J against 
the general hypothesis. This case is then ruled out. 

• If {faipiipi} = 102^20), with \tp) ^ |t/> 2 ), after ap- 
plication of the non-sing ular we obtain the 
coefficient matrix for a generic vector in 2U 

/ aip aipi +13 (3 0\ , 
\ (3 0J [W> 

If aipi + f3 = 0, then there will exist non-null a, f3 
such that a 03^ vector will belong to 2U, against 
the hypothesis. Thus it must be ip\ = 0, i.e. 
\tp) = |0). Under this restriction and after using 
theorem [1] upon the coefficient matrix (jTDJ), all vec- 
tors will be of type W, since r(CW) = 2 for all 
i=l,2, 3, r(M / i) = 2, r(W 2 ) = 1 and the spectrum 
ct(M / 1 ~ 1 W 2 ) is always degenerate. By permutation 
symmetry the cases I^xViV'i) = \<h<p4>2), I^V^), 
with the corresponding restrictions, are also con- 
sidered in this analysis. 



• If \4>i(piipi) ~ \4>2 l P'>p2) , a similar argument leads 
to the coefficient matrix for a generic vector in 2U 
given by 

/ atp + (3 [3 Qaipi\ . . 

\(3 J [ ' 

After application of theorem [TJ it should be clear 
that it is always possible to find non-null a, [3 such 
that (fTTj) corresponds to a GHZ vector in 20, 
against the hypothesis. By permutation symme- 
try, the rest of cases with the 000 generator having 
two factors in common with some of components of 
the W generator is also contained in this analysis. 

• If \<t>itp\il)\)_ = |0iW>), with \ip) ^ |^2,<^2> and 
|?/>) 7^ 1^2,^2)) after a similar argument the coeffi- 
cient matrix of a generic vector in the right singular 
subspace will be 

(atpoipo a(p tpi + (3 a(piip +/3 a(pxipi\ 

^ f3 J [LZ) 

It should be clear that r(CW) = 2 for all i = l,2, 3. 
Now if fiipo + VoV'i 7^ 0, there will always exist 
non-null a, (3 such that r(W / i) = r(M / 2 ) = 1, hence 
a GHZ vector, against the hypothesis. On the con- 
trary, if Lp\i}jQ + ipoipi = Oj then r(Wi) = 2 and 
r (W^ 2 ) = 1, but it is always possible to find non- 
null a, (3 such that ^(PF-f 1 W 2 ) is non-degenerate, 
hence GHZ £ 20. Notice that at most another W 
vector, apart from the generator, can be found. By 
permutation symmetry and by using a x upon each 
qubit, the rest of cases in which the 000 genera- 
tor contains one common factor vector with the W 
generator is also considered in this analysis. 

• Finally, if \(j)iipi?pi) — \<j)ipip), with no common fac- 
tor vector with the W generator, a similar argu- 
ment leads to the coefficient matrix for a generic 
vector in the right singular subspace given by 

/ a</) ip ijj a^oVoV'i + (3 afoipiipo + (3 afoipiipA 
\a(j)itfQipQ + (3 afaipoip! afaip^o a^i^i^iy 

(13) 

A systematic application of theorem Q] makes clear 
that there always exist a, (3 such that this matrix 
corresponds to a GHZ state, against the hypothe- 
sis. 

7. 2U = span{0i*23,0i*23} 

The generic case will be 20 = span{|</>i\I/i), |0 2 vE , 2)}, 
where l^k) denotes a bipartite entangled vector. Several 
cases appear: 
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i. When |0i) and |<fo) arc linearly dependent, then 
2B = \4>) ® span{|*i), |* 2 )}, driving us to the 
already known degenerate classes 0102*, O2GHZ 
and O2IF, respectively, when it adopts the struc- 
tures span}|0*)}, \4>) <g> span{00, 00} and |0) ® 
span{00, *}. In the forthcoming cases, we will then 
impose the linear independence of |0i) and |02). 

ii. If I*!) and 1*2) arc linearly dependent, then 2U = 
C 2 <X> |*), corresponding to the degenerate class 
v E f i2^I / 34, as stated above. Therefore, we will also 
impose the linear independence of 1*1) and | \I/ 2 ) • 

hi. When span}|*i), |* 2 )} = span{00, 00}, then it 
is always possible to hnd linear independent fac- 
tor vectors such that |*j) = a.ji\ipij>) + a^l^^), 
with Oji • cij2 7^ 0. Defining non-singular F^ 
such that FWlaulfa)] = |0), [o 2 i|0a>] = |1>, 
FW{\<p)] = |0), F^m = FW[m = |0), 
F[ 4 1[|V')] = |1) and F^ as in preceding cases, we 
obtain the canonical matrix 



8. 2H = span{0i*23, 02*13} 



1 Ai 
1 A 2 



(14) 



where we have defined A, = ^r-. This matrix cor- 
responds to the canonical state 1 0000) + 1 1100) + 
Ai|0011) + A 2 |llll). with Ai ^ A 2 . 

iv. When span}|*i), |* 2 )} = span}00, *}, then it 
is always possible to find linear independent 
factor vectors such that |*j) = aji\(ftp) + 
cij2 (\<pip) + \<fiip)) , with dji ■ aj2 7^ 0. Defining 
non-smg ular F^ such that F^ [an|0i)] = |0), 
F^[a 21 \<h)} = |1), F®[\ip)] = |0), F^m = \1), 
F [i] [\ip)} = |0), F^[\tp)] = |1) and F™ as in pre- 
ceding cases, we obtain the canonical matrix 



The generic case 
span{|0i*i), 102^2-02) 4 
ties arise: 



is given 

102 V^)}- 



by 2D 

Two possibili- 



i. span{|*i), \1p2ip2)} = span{00,00} = 

span{|i J 92'02), \<f2ip2)}; where the double bar 
also indicates linear independence with re- 
spect to ^2), although possible linear depen- 
dence with 1^2) • Then a generic vector in 2U 

will be given by a|0i) (a\<P2ip2) + &|<^2^2)J + 

P (1 02 (^2^2) + 102^2^2)), with a ■ b 7^ 0. Defining 
non-singular F^ such that F^H^)] = |0), 
F[ 2 ][|02)]-|l),F[ 3 ][|^)]-|O),F[ 3 i[i|^ 2 )] = |l), 
FW[|V 2 )] = |0), FM[|$j)] = |1) and F™ as 
in preceding cases, we obtain the canonical 
vector |0</)00) + |O01V>) + 1 1000) + |1101), with 
the restrictions that \tp) 7^ and \<f>) 7^ |0, 1), 
simultaneously (otherwise we would be recon- 
sidering previous classes). With the use of a x 
upon qubit 4, we are also embracing the case 
in which span} |*i), \(fi2 $2)} = span}00,00} = 
span{|^2' i /'2)j I^V^)}) where now the double_bar 
denotes linear independence with respect to \tp2)- 

ii. span}|*i), \ip2tp2)} = span}00,*} 

span} I ^2) j I ^2^2) + (^2^2)}, with the same 
convention as before for the double bar. Expressing 

|*i) = 01^2^2) + b ^2^2) + Wi^ih and with a 
similar argument to the preceding case, we arrive at 
the canonical state 1 000?/-) + 1 0010) + 1 1000) + 1 1 101) , 
with the restrictions that l^) 7^ and |0) 7^ |0, 1), 
simultaneously (otherwise we would be re- 
considering previous classes). By using <r x 
upon qubit 4, we are also considering the case 
span}|*i), |V>202)} = span}00, *}. 

9. 2n = span}0i*23, 03*12} 



1 Ai Ai 
1 A 2 A 2 



where wc have defined A, 



(15) 



This matrix cor- 



responds to the canonical state |0000) + 1 1100) + 
Ai|0001) + Ai|0010) + A 2 |H01) + A 2 |1110). with 
Ai 7^ A 2 . 



Although this case can also be analyzed with similar 
arguments, we will use permutation symmetry among 
qubits 2 and 3 in the case span}0i*23, 02*13} to arrive 
at the canonical states |O0OO) + |O0<pl) + 1 1000) + |1110) 
and |O0<pO) + |O0O1) + |1000) + |1110). 



10. 2n = span{0 2 *i3, 02*13} 



All possible cases have been considered. Furthermore, 
notice that several states obtained by permutations in- 
volving qubit 1 in preceding cases are enclosed in this 
class. This will also happen in forthcoming classes: some 
permutations involving qubit 1 usually involve a change 
of entanglement class, although the general structure of 
the canonical state shows permutation symmetry [l2| . 



Under permutation symmetry between qubits 2 and 
3 in the case span}0i*23, 01*23}, we rapidly find two 
genuinely entangled canonical states, namely |0000) + 
|1010) + Ai|0101) + A 2 |im) and |0000) + |1010) + 
Ai|0101) + Ai|0001) + A 2 |1011) + A 2 |1110), and one de- 
generate canonical state of type *i3*24 (apart from the 
already considered). 
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11. 2U = span{0 2 * is, 03^12} 

Under permutation symmetry between qubits 2 and 
4 in the case span{0i*23, 03*12} we rapidly find two 
genuinely entangled canonical states, namely 1 000-0) + 
|O010) + 1 1000} + 1 1101) and |O0O0) + |OO10) + |1000) + 
11101). 



12. 2rr = span{0 3 * 12, 03*12} 

Under permutation symmetry between qubits 2 and 
4 in the case span{0i*23, 01*23}, we rapidly find two 
genuinely entangled canonical states, namely |0000) + 
|1001) + Ai|0110) + A 2 |llll) and |0000) + |1001) + 
Ai|0010) + A0O1OO) + A 2 |1011) + A 2 |1101), and one de- 
generate canonical state of type *i4*23 (apart from the 
already considered). 



13. 21J = span{0i* 23 , GHZ} 

The generic case will be of the form 2D = 
span{|0i*i), |02^202) + 1020202)}- Finding non- 
singular transformations as before, we can write 
2D = span{|0*), |000) + |111)}. A generic vector in 
2D will then be of the form a|0*) + /3(|000) + |111>). 
The restrictions on |0) and |*) (thus on |0i) and |*i)) 
will be that no other 0i fc * must belong to 2D. For 
instance, a particular restriction will be that |*) ^ 
span{00,ll}, since otherwise |*) = o|00)+6|ll), a-b 
and then a generic state in 2D will adopt the form 
(cra|0) + /3|0)) |OO) + (a6|0) + 0|1)) |11). It is always pos- 
sible to find non-null a, (3 such that the preceding state is 
a 0i5 , 23 vector different to |0*), i.e. another 0i*23 vec- 
tor: we would be in the span{0i*23, 01*23} case. The 
canonical state will be given by |O0*) + |1000) + 1 1111) . 



14. an = span{02*i3,G#Z} 

By permutation symmetry between qubits 2 and 4, the 
canonical state will be |O* 2 4 03) + |1000) + |HH), where 
the first factor denotes an entangled state between qubits 
2 and 4 and two factor states involving state 1 and 3. 



15. an = span{0 3 *i2, GHZ} 

By permutation symmetry between qubits 2 and 3, the 
canonical state will be |O* 2 3 04) + |1000) + |HH), where 
the first factor denotes an entangled state between qubits 
2 and 3 and two factor states involving state 1 and 4. 



16. an = span{0i* 23 , W} 

By 2D = span{0i*23, W} it is meant that only one de- 
generate tripartite vector of type 01*23 belongs to 2D, all 
the rest being of type W. It is immediate to prove that 
this is impossible, i.e. that there is always possible to find 
a GHZ vector in 2D. The generic case will be given by 

2D = Span! 1 01^100, |02<^202) + 1020202) _+ 102^202)}- 

Defining so that F^\i/j 2 ) = |0), Fl 2 ]|0 2 ) = |1) and 
similarly for k = 3, 4 and using the same F^ as in preced- 
ing cases, the coefficient matrix for a generic [l3[ vector 
|0*) + A (|000) + |111>) in 2D will be 

( 0o*oo 00*01 +A 0o*io + A o *n\ (m) 
^0i*oo + A 0i* O i 0i*io 0i*n0 1 ' 

Now there is a continuous range of values A € C such 
that r(C«) = 2. Focus upon the set of which such that 
detW-2 0= (again a continuous range of them), then it 
is clear that the spectrum of W 2 ~ 1 W\ will be in general 
non-degenerate, driving us to a GHZ vector. Thus we 
are in fact in the case span{0i*23, GHZ}. 

17. an = span{0 2 *i3, W} 

By permutation symmetry the same analysis as in the 
preceding case allows us to conclude that we are in the 
span{02*i3, GHZ} case. 

18. an = span{03*i2, W} 

By permutation symmetry the same analysis as in the 
preceding case allows us to conclude that we are in the 
span{03*i2, GHZ} case. 

19. an = spsm{GHZ, W} 

In this class there exist no degenerate tripartite vec- 
tors, i.e. all of them are either of type GHZ or of type 
W. With the same arguments as before, the canonical 
state is found to be of the form 1 10(^0) + | 1000) + | 0001) + 
|0010) + |0100), where the product vectors |0), |0), \ip), 
|0), |0) and |0) are restricted to produce no degenerate 
state in 2D. 



20. an = span{GHZ, GHZ} 

The notation is meant to describe the situation in 
which only GHZ vectors are contained in 2D. Apart 
for particular instances in which a degenerate tripar- 
tite vector may appear as a linear combination of both 
GHZ generators (which we rule out, since they have 
been already dealt with), we will show that there will 
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always exist at least one W vector in 20. The generic 
case is_ given by span-j^i^i^i) + l^aPaV'a) + 

|02<^2V'2)}- Then, after defining non-singular transfor- 



J 



mations in a similar fashion as in all preceding cases, 
the coefficient matrix for a generic vector in 20 will be 
given by 



r 



(17) 



There is a continuous range of vales A such that 
r(CW) = 2 for all i = 1, 2, 3. Focus upon the set of which 
such that r(Wi) = 2, then there will always exist a value 
A* such that a{W^ W?) is degenerate (A* is the solution 
to the equation A(W r 1 ~ 1 W 2 ) = 0, where A(A) denotes the 
discriminant of the characteristic polynomial of the ma- 
trix A). We are then always in the case span{ GHZ, W}. 



21. W = spanlW, W} 



ticular instances in which a degenerate tripartite vector 
may appear as a linear combination of both W gener- 
ators (which we rule out, since they have been already 
dealt with), we will show that there will always exist 
at least one GHZ vector in 2U. The generic case is 
given by span{|0i<pi'0i) + |0i£it/>i}, |^i^iV'i)> 102^2^2) + 
102^2^2), |</>2</ 5 2'02)}- Then, after defining non-singular 
transformations in a similar fashion as in all preced- 
ing cases, the coefficient matrix for a generic vector in 20 
will be given by 



The notation is meant to describe the situation in 
which only W vectors are contained in 20. Apart for par- 



J 



0iVo'0o+<# , i0o'/ , o+<# , i¥'oV , o+A 0ivoV , i+0iVoi/'i+0i¥'o'/'i ^lViV'o+^i ^lV'o+^iViV'o 0i¥'i'0i+0i | £i>/ , i+0i¥'iV , i 



(18) 



r 



There is a continuous range of vales A such that 
r(C«) = 2 for all i = 1,2,3. Focus upon the set of 
which such that r(M^ 2 ) = 2, then there will always exist 
values A* such that <r(W r 1 _1 W2) is non-degenerate (A* is 
not the solution to the equation A(M / 2 " 1 W r 1 ) = 0, where 
A(A) denotes the discriminant of the characteristic poly- 
nomial of the matrix A) . We are then always in the case 
span{G#Z, W}, too. 



C. Counting entanglement classes 

The preceding classification of states under the SLOCC 
criterion reveals a block structure for the entire state 
space (C 2 ) , with a continuous range of non-equivalent 
vectors inside each block (except for some of them) . Con- 
sidering each block as an entanglement class itself, then 
we have found 18 degenerate and 16 genuine classes, i.e. 
a total of 34 entanglement classes. 



Furthemore, taking into account permutation symme- 
try among the four qubits, indeed we have proven that 
four qubits can be entangled in 8 inequivalent ways. This 
is in apparent contradiction to the result found in [§[, 
where four qubits are claimed to be entangled in 9 pos- 
sible ways. A closer look at this work reveals that one 
of the blocks in that classification contains no genuinely 
entangled vector, namely the class £o 3eI o 3ffl i with canon- 
ical state 1 0000) + |0111). We would like to express our 
conviction that this difference stems from a distinct way 
of counting or considering entanglement classes. We be- 
lieve ours shows a greater agreement with the philosophy 
of the seminal paper by Diir et al [||. Thus we have 
proven the following 

Theorem 2. Let |$) be a genuinely entangled A-partite 
pure state. Then |<f>) belongs to one of the following 8 
genuine entanglement classes and is equivalent to some 
of its canonical vectors: 



The canonical vectors obtained by permutation among qubits 2,3,4 have been omitted in each entry. Restric- 
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Class (2U) 


Canonical States 


Name 


Proposed Notation 


span{000, 000} 


|0000) + |1111) 


GHZ 


2Uooo,ooo 


span{000, fc *} 


|0000) + 1100) + |1111) 
|0000) + |1101) + |1110) 


- 




span{000, GHZ} 


|0<tV^) + 1000) + |1111) 




2Uooo, GHZ 


span{000, W} 


|1000) + |0100) + |0010) + |0001) 


W 


2$ooo,w 


span{O fe *, 0*,*} 


|0000) + 1100) + Ai|0011) + A 2 |llll) 

|0000) + 1100) + Ai|0001) + Ai|0010) + A 2 |1101) + A 2 |1110) 


$4 


2Ho fc *,o fc * 


span{0i*,0j*} 


|O0OO) + lo^iv) + 1000) + |1101) 
\0(f)0tp) + |O01O) + |1000) + |1101) 


- 


2JJo;*,o.,* 


span{O fc <f, GHZ} 


|O0*) + |1000) + |1111) 






span{GHZ, W} 


|0001) + |0010) + |0100) + \l<f><pip) + \l$<prp) 







TABLE I: Genuine entanglement classes for four qubits 



Hons upon the range of the free parameters in each class 
amounts to not producing a precedent structure for the 
right singular subspace 2XJ. 

Names should be given to all classes. A system- 
atic notation for these classes could possibly be 2U$ i ^ , 
where the subindiccs indicate the generators of 2D. In 
this convention GHZ = 2Dooo,ooo, W = 2Uooo,w an d 
$4 = 2Uoi*,Oi*- However a further discrimination refer- 
ring to the particular canonical state must be explicitly 
stated. 

IV. CONCLUDING REMARKS 

We have given the classification of entangled states of 
four qubits under stochastic local operations and classical 
communication using the inductive approach introduced 
in 0] ■ This approach is based on the analysis of the struc- 
ture of the right singular subspace 20 of the coefficient 
matrix of the state in a product basis according to the 
entanglement classes for three qubits. 

Agreeing to consider each structure of 2U as a sin- 
gle entanglement class, we have found 18 degenerate 



and 16 genuine classes (totally 34 classes), where per- 
mutation is explicitly included in the counting. Taking 
into account the permutation among the qubits, 8 gen- 
uine classes can be considered, compiled in table HI As 
expected, in most of the classes a continuous range of 
strictly non-equivalent states is contained, although with 
similar structure. 

This result allows us to predict that there will be at 
most 765 entanglement classes (permutation included) 
for 5-partite systems, 595 at most genuine and 170 at 
most degenerate (cf. Q). No doubt, this classification 
stands up as a formidable task that our inductive method 
may help to elucidate. 
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